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■ Abstract 



We point out two extensions of the relation between matrix models, topological strings 
and A/" = 1 supersymmetric gauge theories. First, we note that by considering double 
scaling limits of unitary matrix models one can obtain large duals of the local Calabi- 
Yau geometries that engineer Af = 2 gauge theories. In particular, a double scaling limit 
of the Gross- Witten one-plaquette lattice model gives the SU (2) Seiberg-Witten solution, 
including its induced gravitational corrections. Secondly, we point out that the effective 
superpotential terms for A/" = 1 ADE quiver gauge theories is similarly computed by large 
N multi-matrix models, that have been considered in the context of ADE minimal models 
on random surfaces. The associated spectral curves are multiple branched covers obtained 
as Virasoro and VF-constraints of the partition function. 
July, 2002 



1. Introduction 



In III we have shown how the effective superpotential in A/" = 1 supersymmetric 
gauge theories, that are obtained by breaking an A/" = 2 super Yang-Mills theory by 
adding a tree-level superpotential for the adjoint scalar $, can be computed by 

a large hermitian matrix models. More precisely, the effective superpotential of the 
N = 1 theory considered as a function of the gluino condensates Si is — apart from the 
universal Si log(5'i/A) terms coming from the pure M = 1 Yang-Mills theory — given exactly 
by a perturbative series that is computed by the planar diagrams of the matrix model 
with potential Furthermore, the contributions of this gauge theory to the induced 

supergravity corrections R^F^^~^ (with R the Riemann curvature and F the graviphoton 
field strength) are similarly computed exactly by the genus (7 > matrix diagrams. 

This gauge theory/matrix model correspondence was a consequence of the large 
dualities of 0,^,^ that relate the computation of holomorphic F-terms in the world- 
volume theories of D-branes to partition functions of topological strings in local Calabi-Yau 
geometries — a relation that was further explored in |[5|,|6|,[7|,^P, p!0| , p!T| . In the simplest case 
these local non-compact Calabi-Yau manifolds take the form 

vv' + y'^ ~W\xf + f{x) =0. 

One finds that in the B-model topological string the tree-level free energy can be computed 
in terms of the periods of the meromorphic differential ydx on the associated Riemann 
surface 

y'^ -W\xf + f{x) = Q. 

As we argued in ^ this curve and the associated special geometry arises naturally 
from the large dynamics of the matrix integral with action But we should stress 

again that the relation with matrix models goes beyond the planar limit. The higher genus 
string partition functions J^g and the related gravitational couplings of the gauge theories 
are exactly computed in the expansion of the matrix models. 

Let us briefly summarize these connections, for more details see [|1| . We start with the 
hermitian matrix integral 



with action 

Sm = — TrM^($), 
9s 
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and W{x) is a polynomial of degree n + 1. The matrix integral can be reduced to an 
integral over the eigenvalues xi, . . . , xat of $ in the potential W{x). In the classical limit 
Qs 0, where one ignores the interactions among the eigenvalues, the equation of motion 
is given by 

y{x)=gsg^ = W'{x)=0. (1.1) 
The associated classical spectral curve is 

/ _ w'{xf = 0, (1.2) 

where x,y can be considered as complex variables. Writing W'{x) = Ylii^ ~ '^i) 

that this singular genus zero planar curve has n double points at the critical points x = at. 

Sometimes it can be helpful to think of the (x,y)-plane as a phase space, with y the 
momentum conjugate to x, as given by the Hamilton- Jacobi equation (|1 . 1|) . Then ( p..2|) has 
an interpretation as the zero-energy level set of the (bosonic part of the) supersymmetric 
quantum mechanics Hamiltonian associated to the superpotential W{x), and S{x) can be 
thought of as the semi-classical WKB action of the associated quantum mechanical ground 
state -^{x) ~ e--^^^^). 

Classically, the eigenvalues will cluster in groups of Ni in the critical points 
where they will form some meta-stable state. The relative number of eigenvalues or filling 
fraction of the critical point we will denote as 

= m/N. 

If Qs is not zero, we have to take into account the Coulomb interaction that results from 
integrating out the angular, off-diagonal components of the matrix $. The equation of 
motion of a single eigenvalue x in the presence of the Dyson gas of eigenvalues xi, . . . ,xn 
is now modified to 

^ 1 

y = W'{x)-2gsYl • 

^ X — Xj 

We will now take the large 't Hooft limit keeping both fj, = QgN and the filling 
fractions Ui fixed. In this case each critical point has its own 't Hooft coupling 

IJ-i = 9sNi = iiVi. 

The collective dynamics of these eigenvalues in the large N limit can be summarized 
geometrically as follows. Each of the n double points x = ai gets resolved into two 
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branch points al^a~ . The resulting branch cuts Ai = [a~,a^] are fiUed by a continuous 
density of eigenvalues that behave as fermions and spread out due to the Pauli exclusion 
principle. This process of splitting up of double points is very analogous to transition 
from the classical to the quantum moduli space in the Seiberg-Witten solution of A/" = 2 
supersymmetric gauge theories [0 — a relation that was explained in [^. The resolution of 
double points is captured by deforming the classical spectral curve ( |1.2|) into the quantum 
curve 

y^-W\xf + lif{x) = Q, (1.3) 

where the quantum deformation f{x) is a polynomial of degree n — 1. The filling fraction Vi 
is related to the size of the branch cut Ai . Roughly, the higher the proportion of eigenvalue 
at the critical point a^, the larger the cut. More precisely, we have 

Hi = liUi = (f y{x)dx 

(This equation can be considered as analogous to the Bohr-Sommerfeld quantization con- 
dition.) Finally the tree-level free energy J-'{fJ.i) can be computed in terms of the dual 
-B-periods by the special geometry relations 



JBi 

where the cycles Bi run from the branch cuts to some cut-off point at infinity. 



2. Unitary matrix models and the Seiberg-Witten solution 

We have reviewed how the effective superpotential in A/" = 1 supersymmetric gauge 
theories obtained by breaking an A/" = 2 super Yang-Mills theory can be computed by a 
large hermitian matrix model. This raises the question whether there exists a matrix 
model that computes directly holomorphic F-terms in the underlying undeformed M = 2 
theory as described by the rigid special geometry of the Seiberg-Witten solution [|12[. As 
explained in |T^] it is indeed possible to extract the geometry of the M = 2 solution from the 
effective superpotential of the M = 1 deformation of A/" = 2 theory and thereby, indirectly, 
from the associated matrix model. But here we will pursue a more direct relation — we 
will show that indeed there is a large matrix model that computes the SW solution 
directly. The matrix model for M = 2 supersymmetric SU{2) gauge theory turns out to be 
a double scaling limit of the most simple unitary matrix model — the so-called one-plaquette 



or Gross- Witten model [|T3 
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2.1. Geometrical engineering of Af = 2 theories 

Our strategy will be the following. The M = 2 theory can be geometrically engineered 
by taking a suitable limit of type IIB string theory on a local Calabi-Yau [^|T3|,|T^,0 . 
This local CY produces directly the Seiberg-Witten curve that encodes the dynamics of 
the Af = 2 gauge theory. More precisely the genus zero topological B-model amplitudes on 
the local CY capture the Seiberg-Witten geometry. The higher genus amplitudes compute 
the contributions of the gauge theory to certain gravitational terms of the form R^F^^~^ 
1^,|19[. We will now engineer a matrix model that is large dual to this local CY 



geometry. In particular the planar limit gives the SW geometry and the corrections 
capture the generation of the corresponding gravitational terms. 

To be specific, let us discuss here the simplest case of pure SU {2) Af = 2 Yang-Mills 
theory. The SW solution is given in terms of the familiar elliptic curve 

w^ = {y^ + uY-K\ (2.1) 

where u is the coordinate on the moduli space, i.e. the vev of the adjoint ^(tr$^), and A 
the gauge theory scale, that we will sometimes set conveniently to A = 1. 

The local CY obtained in the geometric engineering is given by the algebraic variety 

11,1111,113 



vv' -f A^ + -) + 2(y2 + u)=Q 

z 

with z G C*, i.e. an invertible variable. After reducing over the (f , f ')-plane the associated 
Riemann surface is 

A2(z + i) +2(2/2 + ^,) = 0. (2.2) 

Since 2; 7^ we can multiply by z and substituting w = fs?z + j/^ -f w to bring the curve in 
the form 

Furthermore, the reduction of the holomorphic three- form on the CY gives directly the 
SW differential ydz/z. The prepotential J^{u) is then obtained by computing the periods 
of this meromorphic one- form along the A-cycle and S-cycle of the elliptic curve ( |2.1| ). 
Note that there are four branch points at y = ±v— ti±A^. In the classical limit A ^ 
they coalesce pairwise in two double points y = ±^y—u. The moduli space contains two 
singularities at w = ±1 where monopoles respectively dyons become massless. 

Since 2 is a C* variable, it makes sense to write z = e*^ with the variable x periodic 
modulo 2%, and reexpress the original local CY and the resulting curve ( pT2|) as 

A cos X + y + u = 0, 
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with SW differential ydx. This way of writing the equation suggests a relation to a matrix 
model with some suitable potential W{x). In fact, since x is now a periodic variable this 
suggest a unitary matrix model where z = e*^ will get interpreted as the eigenvalue of a 
unitary matrix U. 

2.2. Unitary matrix models 

Unitary matrix models are defined as integrals over the group manifold U{N) of the 

form 

J [ dU-eM--T^rWiU)), (2.3) 

where dU is the Haar measure. As in the hermitian matrix models, one can diagonalize U 
and express everything in integrals over its eigenvalues 

[/~diag(e^"i,...,e^"~), 

where the ai are periodic variables, giving 

^ I KJ I 

Note that such a unitary matrix model can be viewed as a special case of a hermitian 
model by writing U = e*"^ with $ a 'compactified' hermitian matrix, i.e. a matrix with a 
periodic spectrum $ ~ $ + 27r. Such a periodicity is achieved by adding multiples of 2% to 
the eigenvalues a/ of $. This addition of multiple images is the familiar way to compactify 
transverse directions for D-branes or for matrix models in M-theory [^. For example, in 



this way the Vandermonde determinants in the measure become after regularization 

Y[ [ai -cxj + 27rn) = sin(— ^-^-^). 
nez 

The unitary matrix model describes a collection of particles in a potential on 
the unit circle interacting through a Coulomb potential. The equation of motion of the 
unitary model is 

J 

The large solution proceeds exactly as in the uncompactified case. One introduces again 
a resolvent, this time defined as 



, , 1 .X — aj. 



N 
I 
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that satisfies a quadratic loop equation that can be derived in exactly the same way as in 
the hermitian case. In the limit N ^ oo with QgN = fi fixed this loop equation takes the 
following familiar form, when written in terms of the variable y = W'{x) + 2ij.u){x), 

y'^ ~W'{xf + 4fif{x) =0. (2.4) 

The quantum correction /(x) in this unitary case is given by the expression 

/(^) = ^ - ^'(«^)) (2-5) 

I 

2.3. The Gross- Witten model 

Our candidate unitary matrix model will be a much studied one, namely the so- 



called Gross- Witten model [|T3[ with potential VF(a) ~ cos a. This model was originally 
introduced as a lattice discretization of two-dimensional (non-supersymmetric) Yang-Mills 
theory. In such a lattice model to each plaquette with holonomy U around the edge one 
associates the Wilson action (in other contexts known as the Toda potential) 

S(U) = —Tr{U + U-^) = —Tt cos($) 
'29s 9s 

Here $ can be thought of as the lattice approximation to the gauge field strength F^,^, 
and in the limit $ ^ this gives the quadratic Yang-Mills action Tr$^. (The parameter 
e we introduce for convenience. It can of course be absorbed by rescaling Qs = 9ym)- ^ 
general lattice model one integrates over a collection of plaquettes, but in two dimensions 
a single plaquette suffices to compute for instance a Wilson loop action. 
The GW model has two critical points 

W'{a) = -esin(a) = 

at ai = and 02 = tt. Note that for real and positive e (the case relevant for Yang-Mills 
theory) the second point is an unstable critical point. But that issue is irrelevant for the 
holomorphic matrix models that we are considering here. The parameter e can be complex, 
and our eigenvalues e*" are allowed to move off the unit circle into the punctured complex 
plane C*. In fact, following our general philosophy, we will consider the perturbative 
expansion of the matrix integral (|2.3| ) around the saddle point where A'^i eigenvalues are at 
the first critical point ai and N2 = N — Ni are at the second point 02. So we are dealing 
with a two-cut, meta-stable solution to the matrix integral. These two cuts introduce a 
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second parameter, besides the overall 't Hooft coupling n = gsN, namely the relative filling 
fraction 

jy={Ni-N2)/N. 

If we introduce the separate 't Hooft couplings for the two critical points 

fXl = QsNi, fl2 = 9sN2, 

then the difference of these couplings is related to the filling fraction 

n' = fi^- jj.2= gs{Ni -N2) = ^i-v. 

We will be interested in computing the planar limit of the free energy JF as a function of 
the two coupling /ii and or equivalently as a function of the 't Hooft coupling and 
the filling fraction u. Although v takes values in the interval [—1,1] the final result will 
turn out to be a holomorphic function of z/. 

The planar limit can be computed solving the loop equation (|2.4| ). To this end we 
have to compute the quantum correction f[x) defined in (|2.5| ). For our choice of potential 
W{x) = ecosx this becomes the following average over the eigenvalues 



f{x) = - ^e(sinx - sinaj) cot(^^--^) 
I 

= — — : e(cos X + cos a/) (2-6) 



= — e(cosa; + u) 
Here the constant u is defined as the average 

u = — cos ar. 
I 

In the semi-classical approximation gr^ — > we have 

Ai - A2 

u ~ = V. 

A 

since there are Ai eigenvalues at the critical point = and A2 eigenvalues at a/ = tt, 
that contribute respectively +1 and —1 to the average of cosaj. 

Inserting our expression for f{x) into ( ^^.4] ) gives the spectral curve 

— sin^ x + 4|Ue(cos x + = 0. (2-7) 
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We see that the original double points ai = and a2 = tt now split up in four branch 
points , a^. The two branch cuts 

Ai = [aJ",a^], A2 = [a^,aJ] 

describe the condensation of eigenvalues around the two critical points. The discontinuities 
in y give the eigenvalue density Integrating the one-form ydx around the branch cuts 

Ai gives the filling fractions Ni/N . These relations allow one to express the filling fraction, 
or more precisely the relative coupling /j,' = fiu, in terms of the period 

/i' = TT" / yix)dx 
2m J A 

where A is a one-cycle homologous to Ai — A2. This equation describes the exact relation 
between the variables n' and u. Integrating the same differential over the conjugated B- 
cycle encircling the cut [a^^a^'*] computes the variation of the free energy JF under a 
change of the relative number of eigenvalues in the two cuts, 

Note that the GW solution assumes that all eigenvalues center around the stable 
vacuum ai = 0, in which case N2 = and 1/ = 1. In fact we claim that for this vacuum 
we have exactly u = 1, so that the spectral curve is given by 

— e'^ sin^ x 4- 4|Ue(cos x -f 1) = 0. 



Indeed, since one only fills the stable critical point ai = with eigenvalues, only this double 
point will get resolved into two branch points. The second critical point a2 = tt remains 
unresolved. Inserting u = 1 gives the GW solution for the large eigenvalue density that 



can be written as in ||13|| 



~cos(|) J|^-sin2(|). 



2.4- Double scaling limit 

Now we will take a double scaling limit of the GW model to obtain the SW solution 
relevant for M = 2 supersymmetric gauge theory. In this limit we will send N —>■ oo and 
at the same time e ^ and z/ ^ 0, keeping eNi, u/e and Qs fixed — or, equivalently, we will 
send the 't Hooft coupling ^ ^ oo keeping the difi'erence of the two couplings fxi and fX2 



fixed. In this limit the absolute difference in eigenvalues A^i — N2 remains finite, but A'^i 
and N2 become both infinite, and therefore the relative filling fraction u = [Ni — N2)/N 
goes to zero. 

After rescaling y appropriately, the spectral curve reduces in this limit exactly to the 
SW curve (with A = 1) 



Note that the double scaled curve depends only on a single parameter u that at weak 
coupling could be identified with the filling fraction u = [Ni —N2)/N of the matrix model. 
Of course the limit we are taking is at a strong coupling point and the relation between u 
and matrix model modulus 11' is more complicated, as discussed above. 

As we have already mentioned, the prepotential is now computed by the periods of 
the differential ydx = ydz/z along the A and B cycles. This allows us to identify the SW 
periods as 



From the original four branch points , 03 obtained in resolving the two double 
points ai,a2, our double scaling limit takes the branch points a^,a^ to infinity while 
keeping a^,a^ at finite distance. This leaves two homology one-cycles: the S-cycle that 
runs around the cut [a^,a^] and the dual A-cycle that is homologous to Ai — A2. This 
behaviour of the branch points is exactly the behaviour in the double scaling limit one 
takes in the old matrix models [pT| , p^ , as we also noted in [|l|. For example, the (2,3) 
critical point of the one-matrix model was obtained starting from a curve of the form 
y"^ = + . . . and the double scaling limit got rid of the all monomials with power more 
than x^ giving an equation of the form = x'^ + . . . (Note that these branch points in the 



/i' = gs{Ni - N2) = niy 



y'^ -\- cos X -\- u = Q. 




(2.8) 
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x-plane should not be confused with the branch points in the y-plane that are relevant for 
the SW solution.) 

The double scaling is very analogous to the limit that was used in the A-model topo- 
logical string in [^. In fact, when the scaling limit is embedded in type II string theory, 
the resulting CY geometry based on (|2.7| ) will have RR flux through the compact cycles 
Ai and A2. It is crucial that the remaining compact cycle B does not carry any Ramond 
flux. We are thus engineering a large dual of a geometry without fluxes. 

The GW model has a famous third order phase transition at (in our convention) 
n/e = 2. This signals a transition of the eigenvalue distribution in which the single cut 
changes topology and starts to cover the whole unit circle. After the GW phase transition 
the eigenvalue distribution is given by p{a) ~ cosa + p/2. Geometrically speaking, in that 
phase all four branch points are on top of each other. 

This phase transition is however not relevant in our model. First of all, we are studying 
a more general question by considering a stationary phase approximation around a meta- 
stable state with two clusters of eigenvalues and consequently have to work with a two- 
dimensional phase diagram (/x, u). As we argued the GW solution puts = 1 and that is 
very far away from our double scaling limit in which u tends to zero. Indeed in our limit 
the number of eigenvalues in the two cuts is roughly equal. Secondly, we are dealing with 
an holomorphic object, and holomorphy excludes any phase transitions, one can just go 
around the singularity. The GW phase transition is just a (very special) real slice of our 
complex phase diagram. 

Finally it would be interesting to connect this approach to the beautiful semi-classical 



computation of the SW solution and its gravitational counterparts in That compu- 
tation was inspired by matrix integrals appearing in D-branes formulas. 



2. 5. Generalizations for general groups 

It is not difficult to guess how the SW solution for gauge group SU (n) can be engi- 
neered. In this case the curve associated to the local CY is of the form 

cosx + Pn{y) = 0, 

with Pn{y) a polynomial of degree n. More generally we can consider a chain of U{ni) gauge 
theories with bifundamental matter, for which the corresponding curve has been obtained 
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from the M5 brane viewpoint in and from the viewpoint of geometric engineering in 
1% . This wiU give rise to a curve of the form 



i^(e^^?/) = 

where F is a polynomial in e^*^ and y. In particular if we consider the rank of all the 
gauge groups to be equal to n, then F is a polynomial in y of degree n. Moreover the 
difference in power of e*^ between the highest and lowest powers is the number of U{n) 
gauge groups plus 1. As we will discuss in greater detail in the next section in the context 
of hermitian multi-matrix models such curves are typically produced by a multi-matrix 
model consisting of n — 1 matrices. The choice of the coefficients in F will be related to the 
choice of the action and some suitable double scaling limit, as we studied in the context of 
SU{2) gauge theory here. 



2.6. Connections with A-model topological strings 

It is natural to ask if there is a connection with A-model topological strings, and in 
particular for A-models on local toric CY. As was demonstrated in |26,27| the B-model 
mirrors are given by 

W + F(e*",e^^) = 

for some F. This is analogous to an infinite matrix model version of the unitary matrix 
models, as follows from our discussion above. In this case the A-model has a gauge theory 
dual involving certain correlations functions of Chern-Simons theory It would be 

interesting to connect these matrix models directly with the Chern-Simons gauge theory 
computation, thus completing the circle of ideas. There are some hints that this idea 
indeed works 



3. Quiver matrix models 

We will now turn to a related generalization of where we will connect superpotential 
computations of quiver gauge theories to multi-matrix models. 
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3.1. Af = 1 quiver gauge theories and topological strings 

We will restrict our discussion here to the ADE quivers, in particular the case, 
although one can also include the affine quivers based on the extended Dynkin diagrams 
ADE. Let r denote the rank of the quiver G and consider a partition 

N = N^ + ... + Nr. 

In the associated Af = 2 quiver gauge theory we assign to each of the r vertices Vi of 
the Dynkin diagram of G a U{Ni) gauge field and to links connecting vertices Vi and 
Vj we associate bifundamentals Qij transforming in the representation (Ni^Nj) with a 
hermiticity condition qIj = Qji. 

For such a gauge theory we can write a general tree-level superpotential 

= 5^s,,TrQ,,$,g,, + 5^TrVr,($0, (3-1) 

with Sij = —Sji = 1 (for some ordering i < j), if the vertices Vi and Vj are linked in the 
Dynkin diagram (we will write this relation also as {i,j)), and Sij = otherwise. Here the 
first term is the standard superpotential of the M = 2 theory with bifundamental matter. 
The additional potentials Wii^i) are introduced to break the supersymmetry down to 
U=l. 

Within type II string theory these quiver gauge theories are obtained by wrapping 
D5-branes over a particular CY geometry that is a fibration of the corresponding ADE 
singularity over the complex plane 0. This geometry contains r intersecting P^'s. Ac- 
cording to 0,1^ in the large limit the geometry undergoes a transition to a deformed 
geometry where these P^'s are blown down and a number of S'^^s with RR flux are "blown 
up." The corresponding smooth CY geometry gives a dual description of the gauge theory 
system. 

In the context of B-model topological strings, the deformed CY geometry is dual to a 
two-dimensional large gauge system, obtained from a collection of B-branes wrapped on 
the intersecting P^'s. The world- volume theory of these branes consists of open topological 
strings. So each P^ gives rise to a two-dimensional field theory with Lagrangian 



S{^) = - [ TtUIDa^"! + W^.($°V). (3.2) 

gs jpi ^ ^ 
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where u is some volume form on P^, and $^ and $| are adjoint fields of respectively spin 
and spin 1 coupled to an U{Ni) holomorphic gauge field. Here we included the effect 
of the superpotential Wii^o). The open topological strings connecting different P^'s give 
as physical fields the bifundamentals Qij. Since the different P^'s intersect in points, the 
action of these bifundamental scalar fields localizes to the intersection point x and is given 
by 

(Compare the similar computation for the coupling of open topological strings connecting 
Lagrangians A-branes intersection along one- dimensional curves in pi[|.) 

As in |jl|] one can see that in the end this two-dimensional topological field theory can 
be completely reduced to the zero modes of the fields $°(x) = $i and Qij. Thereby the 
path-integral reduces to the "quiver matrix integral" 

Z= J]dQ,,exp(--TrVr($,Q)). (3.3) 

where $i (with z = 1, . . . , r) is an A^i x A^j hermitian matrix and for every linked indices 
the variable Qij is a A^ x Nj rectangular matrix, satisfying Q\- = Qji. 
The generalization of the conjecture in will now identify the free energy of the 
large A quiver matrix model, for given filling fractions of the saddle points, with the 
closed topological string partition function in the corresponding deformed CY geometry, 
and this in turn with the effective superpotential of the quiver gauge theory. 

3.2. Saddle points and dual CY geometry 

The saddle points of the quiver superpotential have been discussed extensively in 
following the mathematical literature. The eigenvalues x of the adjoints $i have to 
satisfy a series of equations: one for every positive roots ak of G. If that root is expressed 
in the simple roots as 

i 

hen the associated condition reads 

E^lVr/(x)=0. (3.4) 

i 
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The saddle points can be labeled as Xa,k with a = 1, . . . ,d and d the maximal degree 



occurring in (|3.4| ). If such a critical point appears with multiplicity Na^k then the total 
number of eigenvalues of the matrix $i in this saddle point is Na^k ' A general saddle 
point is therefore parametrized by the filling fractions z^a,/c = Na^k/N. We will consider the 
matrix integral in the limit where both the A''^ ^ and tend to infinity keeping the filling 
fractions and the 't Hooft coupling finite. 

In the case of an Ar quiver there is a more straightforward description of the saddle 
points . Introduce the r + 1 potentials 

i 

to{x)=0, ti{x) ='^Wj{x), i = l,...,r. 

The the associated rational planar curves 

y - ti(x) = 
intersect in various double points given by 

j 

tj{x)-U{x)= J2 Wl{x)=0. 

k=i+l 

The saddle points of the Ar quiver matrix potential correspond exactly to these double 
points. 

In this Aj. case the original singular CY geometry is given by 



r 

UV 

1=0 



+ ll{y~Uix))^o 



which after reduction over u, v gives precisely this collection of nodal curves. After the 
deformation the corresponding smooth Riemann surface is given by 

r 

l[{y-U{x))+f{x,y)=0. (3.5) 

i=0 

for a suitable normalizable quantum deformation f{x,y). Again every double point gets 
resolved into two branch points. The resulting quantum curve is now an r + 1 fold cover 
of the x-plane. By moving around in the x-plane these sheets will be exchanged through 
Weyl reflections acting on the parameters ti. 
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The analogues of the meromorphic one-form are constructed by reducing the holo- 
morphic three-form of the local CY over various cycles and for Aj. have the following 
description [^^. Write the curve (|3.5|) in the factorized form 

r 

l[{y-a,{x))=0. 

i=0 

Then we have a basis of r — 1 meromorphic one- forms r]i, . . . , rj^-i that is in one-to-one 
correspondence with a basis of positive roots for Ar-i given by 

rji = (ai_|_i(a;) — ai{x))dx. (3-6) 

In the undeformed case, with ai{x) = ti{x), this gives rji = W-{x)dx = dWi. Different 
choices of positive roots correspond to Weyl reflections, which are in fact the generalized 
Seiberg dualities of gauge theories (which should also have some direct interpretation as 
dualities of the matrix models). Note that in the case of Ai we have 

(y-a)(2/ + a) = 0, a'^W'ixf + fix) 

and this gives the usual deflnition 77 = ±ydx. 

We will argue that this process of smoothing out the singular curve including the 
meromorphic differentials is exactly described by the large dynamics of the quiver matrix 
integral. 

3.3. ADE matrix models 

Quite remarkably it turns out that the quiver matrix integrals (|3.3|) (up to some minor 
details) have already been studied in the context of the "old matrix models." They have 
been used to describe the coupling of ADE conformal minimal models to two-dimensional 
gravity by Kostov [|3^, see also the reviews ||3^,0, and they have naturally emerged in 



the study of matrix models and integrable systems in the work of the ITEP group ||35|| . 
We will follow closely these works in presenting the main results, leaving the details to the 
literature. 

First of all, one can immediately integrate out the bifundamental flelds Qij in the 
quiver matrix integral to give an effective interaction between the adjoint flelds $i and $j 

det($i (g) 1 - 1 (g) 
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After expressing everything in terms of the eigenvalues of the remaining hermitian matrices 

$i ~ diag(Ai,i, . . . , Xi,Ni) 
the quiver matrix integral reduces to 

, 2 



z 



n (v-M' 

i<j,I,J 



Note that if we introduce the Cartan matrix of G 

Cij ^^ij I "^ij I ' 1 

then the quiver eigenvalue measure, that generalizes the usual "fermionic" Vandermonde 
determinants of the one matrix model 



, 2 
KJ 



can now be written as 

The resemblance to a correlation function of vertex operators is not accidental — it was in 



fact the main motivation to study these kind of matrix models in ||3^ since it allows one 
to express the partition function as a particular state in a two-dimensional chiral CFT, to 
be more precise the level one realization of the corresponding ADE current algebra. 

In the large N limit this many-flavor Dyson gas of eigenvalues will spread out in cuts 
around the saddle points and will form a continuum of eigenvalues described by a series of 
densities functions 



I 

The solution of the model proceeds again through the resolvents or loop operators of the 
matrices $j 

1 „ . 1 X 1 1 



ojAx) = — Tr( ) = — > . 

The jump of uji{x) across a branch cut measure the eigenvalue density Pi{x). 
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In fact, it is natural to work with a closely related object — the derivative of the matrix 
model action 5" with respect to an eigenvalue of type i evaluated at a general position x 
in the complex plane (away from the cuts) as it interacts with all the other eigenvalues, 



y,{x) = Qs d,S = W^x) - 2i,C'^Uj{x). (3.7) 

As we already mentioned in for a multi-matrix model we want to identify the one-forms 
rji (|3.6|) coming from the local CY geometry with the expressions yi{x)dx, up to a possible 
change of basis. The most powerful and general technique that can be used to relate the 
CY geometry to the matrix model are the loop equations. 

3.4- Loop equations and collective fields 

Before we discuss the loop equations of the multi-matrix models, let us first rewrite 
the solution of the one-matrix model as used in |1|] in a more suggestive form, that is 
actually a standard technique in matrix model technology. Here we found among others 



the reviews [^,^^ very helpful. 

The resolvent uj{x) has a natural interpretation as a loop operator. More precisely, 
the inverse Laplace transform 



f '^e'-'uj{x)=TT{e"^) 



is the zero-dimensional analogue of the Wilson loop. The non-linear all-genus loop equation 
is usually written in terms of uj{x) as 

'-{u{z)) = ^i{u:{xY) (3.8) 



Q 27vi X - 

where ( ■ ■ ■ ) indicates an expectation value within the matrix integral. The contour C 
encircles all the cuts but not the point x. This equation is supplemented with the boundary 
condition ( ^(a;) ) ~ l/x at infinity. The loop equation acts as a Schwinger-Dyson equation 
of the matrix model. It gives a recursive relation to solve for the loop operator and the 
free energy. In the planar limit we have large factorization ( t<;(a;)^ ) = ( c<;(a;) )^ and the 
loop equation becomes algebraic. 

Loop operators are closely connected to collective fields. By integrating out the angu- 
lar variables the individual eigenvalues start to behave as fermions, and the collective field 
is essentially constructed by bosonization of these fermion fields. In |^ we have already 
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speculated that this coUective field should be identified with the Kodaira- Spencer field ^8 
describing the closed strings moving on the local CY geometry. 

For a single matrix model the collective field is defined as the chiral two-dimensional 
scalar field 

i^[x) = W{x) - 2gs ^ log(x - A/) 

It clearly satisfies dip{x) = y{x) with 

y = W'{x) - 2g, V = W'{x) - 2fiu;{x) 

^ x-Xi 

So in view of ( |3.7D we can identify the function ip{x) with the action S{x) of a single 
eigenvalue as a function of its position x in the complex plane in the presence of the gas 
of other eigenvalues Ai, . . . , Xn- The function ^{x) is multi-valued in the x-plane. It has 
branch cuts around which it changes sign. It is therefore only properly defined on the 
double cover 

2/2 - I^'(x)2 + /(x) = 0. (3.9) 

On this Riemann surface ip{x) has quantized periods around the A-cycles, given by the 
filling numbers /i^ = QgNi. Since it is a chiral field the periods around the dual S-cycles 
are not independent and expressed by the special geometry relations as dT ld\Xi. 
Note that if we work with a general, not necessarily polynomial, superpotential 



W{x) = Y,t 



n>0 



then the expectation value of the field dip{x) inserted in the matrix integral can be repre- 
sented by a linear differential operator in the couplings tn acting on the partition function. 
For example, 

(a^(x)) = (J] nt^x^-^ - 2gl ^' 



-1 d 

L 

n>0 n>0 



and similarly for multi-point functions. (Here we used that the derivative d/dtn brings 
down a factor — Tr$".) 

With this notation there is an elegant way to write the loop equations. Introduce the 
holomorphic stress-tensor 

T{x) = {d^f = Y,Lr,x-^-^ 



n 
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Then the aU-genus loop equation (|3.8|) of the one-matrix model can be rewritten in the 
suggestive form 

" dz 1 

'T(^)) = 0, (3.10) 



i 



}(^ 2m X — z 

That is, the expectation value {T{x)) has no singular terms if a; ^ 0. Therefore we can 
also express (|3.10|) equivalently as the Virasoro constraints [^^| 



LnZ = 0, n > —1. 



The derivation of the constraints in the matrix model is completely standard — it simply ex- 
presses the Ward identities following from the invariance under infinitesimal reparametriza- 
tion of $ — s> $ + e^"^"*"^ of the matrix variable $. 

In the planar limit we can substitute the classical values for dip{x) = y m. T[x) and 
then equation ( p.lO| ) is a consequence of (|3.9| ) that can now be written as 

T{x) = W'{xf-f{x), 

which shows that T{x) is indeed regular (even polynomial) at x = 0. 

3.5. Quiver theories and W -constraints 



The large N solution of the quiver matrix integral ( |3.3|) now proceeds along similar 
lines [ p2| , |35| . One introduces r scalar fields ^pi{x) through the one- forms (|3.7| ) as 



yi{x)dx = d(pi{x) 

One can then show that the multi- valued fields ipi{x) are actually the values of one single- 
valued field ip{x) (essentially the full matrix model action) on a r -f 1 branched cover of 
the complex x-plane. This branched cover is the spectral curve associated to the quiver 
matrix integral, and turns out to be given by ( |3.5| ) in the A^. case. 

The general derivation of the curve proceeds through generalized loop equations. 
For these multi-matrix model we do not only have the Virasoro constraints, expressing 
reparametrization invariance in the matrix variables $i. There are also higher order re- 
lations |]3^,^. The full set of loop equations are obtained by showing that the partition 



function Z satisfies a set of W^-constraints, labeled by the Casimirs of the corresponding 
ADE Lie algebra, which contains the Virasoro constraints. These constraints take the 
form 

Q 2m X — z 
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where W^^\x) is a spin s current in the VF-algebra. When expressed in modes these 
equations take the form 

Wi'^ -Z = 0, n>l-s. 

In the case of there is leading spin r + 1 current that with a suitable basis of vectors 
(fo, . . . ,Vr can be written as 

r 

W^^'+^Hx) ^ l[{v, ■ d^) + . . . 

1=0 

We claim that in the planar limit this loop equation translates directly into the curve ( |3.5|) . 

To be completely explicit let us give some more detail for the simplest case of A2. Here 
we have two matrices $1, $2 with potentials Wi{^i) and VF2($2)- The classical singular 
curve is after a shift in y given by 

{y-h{x)){y-t2{x)){y-ts{x)) =0 

with 

t^ = -(^2W[ + W^)/3, t2 = {W[-W^)/3, ts = {W[ + 2W^)/3, 
all polynomials in x. To find the quantum curve we introduce the resolvents 

Wlix)=Y] \ . W2{x)=Y] , 

and the one-forms yi{x)dx 

yi = W[ - ^{2ui - U2), 2/2 = - n{2u2 -ui). 
We now claim that the quantum curve is given by 

{y - ai(x)) [y - 02(0;)) [y - a3(x)) = 0, 
where the functions ai{x) are no longer polynomials, but instead are defined as 

ai=ti-f|UWi, a2 = ^2 - A*(i^i - 1^2), as = ts - IJ.UJ2. 
With this choice we have, as claimed before, 

a2-ai= yi, as - 02 = 2/2- 



20 



Now after some algebra, expanding out terms like u!i{x)^, one verifies that indeed 

{y - ai{x)){y - a2{x)){y - asix)) = 

{y - h (x)) {y - t2 (x)) {y - hix)) + f{x) y + g{x) = 

with f{x) and g{x) polynomials. 
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